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The shear viscosity r\ of a quantum liquid in the vicinity of T\ is examined. In liquid helium 
4 above T\ (T\ < T < 3.7 K), under a strong effect of Bose statistics, the coherent many-body 
wave function grows to an intermediate size between a macroscopic level and a microscopic 
one. These wave functions are qualitatively different from thermal fluctuation, and manifest 
themselves in the gradual decrease in shear viscosity above T\. To formulate this phenomenon, 
we combine the correlation function with fluid dynamics. Applying the Kramers-Kronig relation 
to the generalized Poiseuille's formula for capillary flow, we perform a perturbation calculation 
of the reciprocal 1/r] with respect to the particle interaction, and examine how the growth of 
coherent wave functions gradually decreases shear viscosity. Comparing with the experimentally 
determined r/(T), p B (T)/p of such a mesoscopic condensate is estimated to reach 10 -5 just above 
T\. We examine the effect of condensate size on the stability of such a superflow, and touch 
upon the superflow in porous media. 

KEYWORDS: Superfluidity, Capillary flow, Liquid helium 4, Shear viscosity, Mesoscopic condensate 
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1. Introduction 

Superfluidity was first discovered in a macroscopic fric- 
tionlcss flow of liquid helium 4 through a capillary or a 
narrow slit. 1 ' Since then, the range of studies of super- 
fluidity has expanded. Since the late 1970's, when people 
had an impression that our understanding of superflow in 
bulk liquid helium 4 was essentially completed, 2 ) exper- 
imenters have turned their interest to superflow in more 
exotic structures such as porous media. 3 ' These materi- 
als allow us to study new aspects of superfluidity, so far 
inaccessible in the bulk liquid helium 4. As is often in 
physics, however, the original capillary flow experiment 
that first demonstrated superfluidity is not the most con- 
ceptually clear-cut demonstration of it. Capillary flow 
is well known but not completely understood even at 
present. 4 ' One can think of two reasons: (1) Capillary 
flow is a nonequilibrium phenomenon 5 ' , 6 ' which exhibits 
subtle features that do not exist in the equilibrium man- 
ifestation of superfluidity such as nonclassical rotational 
properties. (2) In supcrfluid helium 4, many features as- 
sociated with Bose statistics are masked by the strongly 
interacting nature of the liquid. The formulation of the 
viscosity of a classical liquid, especially from the struc- 
tural viewpoint, 7 ' has been a difficult problem, because 
complex motions inherent in the liquid do not allow us 
to make a simple microscopic formulation. 

In this paper, the capillary flow above T\ is revisited. 
The two-fluid model, which is the basis for understanding 
liquid helium 4, separates the system into the normal and 
superfluid parts from the beginning, and normally as- 
sumes that the latter abruptly appears at T\ . The success 
of the two-fluid model instilled into us the notion that su- 
perflow appears at T\ in a mathematically discontinuous 
manner, and that all anomalous properties above T\ arc 
due to thermal fluctuation. (In fact, in liquid helium 4 at 
\T/T\ — 1| < 10~ 2 , the sharp increase in specific heat or 
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Fig. 1. The kinematic viscosity v(T) (a solid curve) of a bulk 
liquid helium 4. 



the softening of sound propagation arises from thermal 
fluctuation.) Figure 1 shows the temperature dependence 
of the kinematic viscosity v(T) — rj(T)/p of the capillary 
flow of liquid helium 4 s ' 9 '. 10 ' In Fig.l, when cooling the 
system, the kinematic viscosity v(T) does not abruptly 
drops to zero at T\ like a dotted line, but it gradually 
decreases in T\ < T < 3.7 K, and finally drops to zero 
at T\. Similarly to the anomalies at \T/T X - 1| < I0~ 2 , 
this gradual decrease in v(T) in T x < T < 3.7K has 
been tacitly attributed to thermal fluctuation (except in 
ref.8). Actually, behind the two- fluid model, there ex- 
ists a more basic assumption in statistical physics: the 
infinite-volume limit. To clearly define the phase transi- 
tion, this limit eliminates the intermediate-sized order, in 
which "intermediate size" is equivalent to "microscopic 
one" . For the mechanical properties, however, in con- 
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trast to the thermodynamic ones, one cannot ignore the 
boundary condition of the system. In particular, in this 
system, since the size of the coherent wave function and 
that of the capillary are not so different at low tempera- 
tures, we cannot naively use the V — > oo limit. Rather, 
the magnitude of phenomena masked by this limit should 
be estimated by experiments. (In this sense, the low- 
temperature phenomenon with quantum coherence is es- 
sentially a mesoscopic one.) 

Physically, the interpretation of v(T) in T\ < T < 
3.7K on the basis of thermal fluctuation is questionable 
for the following reasons. (1) The range of temperature 
of 1.5 K from T\ to 3.7 K is too large for the thermal fluc- 
tuation. 11 ) For the fluctuation of temperature AT, one 
knows the formula {(AT) 2 ) = k B T 2 /C v , in which C v is 
the heat capacity of a fluctuating small region. As men- 
tioned in ref.8, the Cy in the case of AT = 1.5K is about 
2.8 x lCr 23 J/AT, which implies that the fluctuating re- 
gion would have to be as small as one atom diameter. (2) 
Thermal fluctuation gives rise to the short-lived and ran- 
domly oriented wave function obeying Bose statistics. In 
contrast, for the one-directional flow to show very small 
shear viscosity, 12 ' a long-lived and coherent translational 
motion of particles along a specific direction is necessary, 
which is qualitatively different from thermal fluctuation. 
(In this sense, the fluctuation-dissipation theorem is not 
applicable to the nondissipative response such as the de- 
crease in viscosity. 13 )) Rather, the gradual fall of v(T) in 
T\ < T < 3.7 K is likely to represent the orderly redis- 
tribution of particle momenta due to Bose statistics , an 
advance sign of the state below T\.^ 

In the Bose system just above T\, particles experience 
a strong effect of Bose statistics, and the coherent wave 
function grows to a large but not yet macroscopic size. 
Since these wave functions are negligible at V — > oo, 
they do not affect the definition of phase transition at 
T\. However, in contrast to thermal fluctuation, these 
intermediate-sized coherent wave functions have a pos- 
sibility of leading to a coherent translational motion of 
particles within a mesoscopic distance. 14 - 1 At T\ < T < 
3.7 K, there is no macroscopic condensate connecting the 
two ends of the capillary, hence no macroscopic friction- 
less transport occurs. Rather, the intermediate-sized co- 
herent wave function is likely to reduce the v(T) of cap- 
illary flow. We will call such a wave function mesoscopic 
condensate not in the thermodynamic sense but in the 
mechanical sense. Figure. 1 must include valuable infor- 
mation on the mesoscopic condensate. Our problem is to 
estimate it by experiment and compare it with the micro- 
scopic model. For this purpose, we must take a somewhat 
different approach to supcrflow. 

(1) The shear viscosity of a quantum liquid has 
been studied by applying the kinetic theory of gases to 
phonons and rotons. 15 ) However, for our purpose, this 
method is not applicable. Well below T\, various exci- 
tations of liquid helium 4 are strictly suppressed except 
for phonons and rotons. Hence, they are normally re- 
garded as a weakly interacting dilute Bose gas, although 
the excitation in liquid. For v(T) above T\, however, the 
dilute-gas picture has no basis, because the basis for this 
picture, the macroscopic condensate, has not yet devel- 
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Fig. 2. The change of the conductivity spectrum <t(lo) from (a) 
in a classical fluid to (b) in a supcrfluid. 



oped. Rather, the strongly interacting excitations in the 
liquid determine the v above T\. The effect of Bose statis- 
tics on such excitations is worth studying. 

In §2, we will develop an alternative method combin- 
ing the correlation function with fluid dynamics. 16 ) The 
distribution of flow velocity in a capillary is described by 
Poiseuille's formula 

d 2 r 2 AP 

where d and L are the radius and length of the capil- 
lary, respectively, and AP is the pressure between the 
two ends of the capillary. The simplest method of mea- 
suring n is to set the capillary to stand vertically at grav- 
ity g. The level h(t) of the liquid decreases as h(t) = 
h(0) exp(— at), where a includes n as pgd 4 /(nL) (see Ap- 
pendix. A). Hence, AP(t) also decreases as AP(t) = 
AP(0) cxp(-crf) in eq.(l). In T A < T < 3.7K, the curve 
of h(t) shows an anomalously small shear viscosity. Al- 
though it shows no sign of frictionless transport, the fre- 
quency spectrum of the quantities in cq.(l) reflects the 
mesoscopic dynamics in the flow as well as the macro- 
scopic one, and therefore, mesoscopic supcrflow must ex- 
hibit a characteristic spectrum in it. On the axis of the 
capillary (r = 0), multiplying both sides of eq.(l) by the 
density p, we define the massflow density j (t) = pv x (0, t) , 
and the conductivity a = p/(4r/). Under an oscillating 
P(uj)exp(iu)t), the spectrum j(u) and er(cj) (Fig. 2(a)) 
satisfy 



j(uj) = -a(u)d 



*AP(w) 



(2) 



In the capillary flow below T\ or in the vicinity of T\, 
the superfluid and normal-fluid parts flow without any 
transfer of momentum to each other. Accordingly, o~(u>) 
splits into a sharp peak A5(w) at u — and a continuous 
spectrum a n (uj) as in Fig. 2(b); 



j(u>) = - [a n (uj) +A5(u)]d 



2 AP( W ) 



(3) 



In §2, eqs.(2) and (3) will be embedded into the gen- 
eral linear-response relation including both the dissipa- 
tive and nondissipative (dispersive) responses. Applying 
the Kramers-Kronig relation to them, we will derive a 
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formula of u{T) = 1/[4<t(T)] (eq.(ll)) in the vicinity of 
T\, and quantitatively estimate the effect of Bose statis- 
tics in Fig. 1. In §2.3, the role of Bose statistics in the 
fall of shear viscosity will be discussed using the Maxwell 
relation. 

(2) In the microscopic theory of shear viscosity 77, one 
normally starts from the two-time correlation function 
of the tensor J xy (t) = — J2 i (Pi,xPi, v / , m), and perform a 
perturbation calculation of (J xy (0)J xy (t)) with respect 
to the particle interaction 5 ). 6 ) This method has the fol- 
lowing difficulties in solving our problem. In the weak- 
coupling system (a gas, or a simple liquid such as liquid 
helium 4), the particle interaction U normally enhances 
the relaxation of particles to local equilibrium positions 
in the flow, thereby decreasing 77. 17 ) If we would try to 
formulate this property using the Kubo formula for 77, 
we must derive the decrease in 77 from the increase in U. 
Since 77 and U changes in opposite directions, there must 
be a delicate cancellation of higher-order terms in the 
perturbation expansion of r\. In contrast, when we start 
from Poiseuille's formula, in which 77 appears in the de- 
nominator of the linear- response coefficient as in eq.(l), 
we apply the Kubo formula not to 77 but to the reciprocal 
I/77 in eq.(2). In the perturbation expansion of I/77 with 
respect to U, the increase in U naturally leads to the in- 
crease in I/77, and therefore to the decrease in 77 . In this 
case, one need not expect the cancellation, and the effect 
of U on 77 is simply built in from the beginning. 

(3) We will focus on the redistribution of particle mo- 
menta in the vicinity of T\. In §3, following ref.14, we as- 
sume the Bose system without the macroscopic conden- 
sate as a non-perturbed state, and perform a perturba- 
tion calculation of I/77. By taking characteristic diagrams 
reflecting Bose statistics in the expansion, we will exam- 
ine how the formation of a larger coherent wave function 
gradually decreases shear viscosity. Specifically, combin- 
ing in §3 the phenomcnological relation eq.(3) with the 
microscopic calculation, we will derive the formula of 
er(cj) (eq.(28)) illustrated in Fig. 2(b) , and derive the 
suppression of v. 

(4) Compared with macroscopic condensate, the meso- 
scopic condensate shows different responses in the oscil- 
lation experiments. After defining the damping angular 
frequency (eq.(24)), we will examine the stability of su- 
perflow to the oscillation, and estimate the effect of con- 
densate size on the stability (eq.(31)). From this view- 
point, we will briefly discuss the response of the superflow 
in porous media to torsional oscillation or ultrasound in 
§4- 

Lastly, we will discuss these results in a wider context 
in §5. 

2. Shear Viscosity in Capillary Flow 

2. 1 Formalism 

In classical fluids, the capillary flow is a typical dissi- 
pative process. To describe both dissipative and nondissi- 
pative processes, let us generalize the conductivity spec- 
trum er(c<j) in eq.(2) to the complex number 

= - [a, ( W ) + 10 2 ( W )] d 2 . (4) 



The spectrum of AP(t) = AP(0) cxp(— at) is given by 
AP(cj) = AP(0)tt- 1 (cj 2 + a 2 )- 1 / 2 , which has a peak 
at uj = 0. We define the half width coj of the peak as 
ujf = \/3a. In eq.(4), <7i(w) must satisfy the following 
sum rule 6 ) including the conserved quantity f(d) (see 
Appendix. B) 

- / ai(w)dw = f(d). (5) 
* Jo 

The average conductivity a is defined as a = 
^f 1 So' <7 i( w ) £ ^ w i hence giving the shear viscosity 77 = 
p/{4a). 

To incorporate eq.(4) into the linear- response theory, 
instead of AP(u), we will use the velocity 7j (w), which 
satisfies the equation of motion pdvo(t)/dt = —AP(t)/L. 
Since vq(u>) and j(ui) constitute the perturbation energy 
/ j(u))v (uj)duj, we can regard v (uj) as the external field 
in eq.(4). We rewrite eq.(4) as 

j(w) = p[-ua 2 (w) + iuja 1 (uj)]d 2 v (uj) 7 (6) 

and regard p [— uja 2 (uj) + iuja\{uj)] d 2 as the coarse- 
grained form of generalized susceptibility. In this for- 
mula, the real part —ua 2 (uj) corresponds to the non- 
dissipative process, and the imaginary part coa\(uj) to 
the dissipative one. In liquid, the former corresponds to 
the Couette flow of viscous liquid in a rotating bucket. 
When the whole fluid is in the uniform rotation, no in- 
ternal friction occurs, 18 ) and it rotates like a rigid body 
as Vd(r) — fi x r (O is the angular velocity). The latter 
corresponds to the capillary flow. They are connected to 
each other by the Kramers-Kronig relation 19 ) as 

2 [°° ua 2 (u) 

Using the response of the rotating liquid as the real 
part coa 2 (u>), we can obtain the conductivity of the cap- 
illary flow ai(co) by eq.(7). u;a 2 (uj) is derived from the 
susceptibility XuviliW) (the real part of the general- 
ized one). 20 ) In general, the X^il^) 1S decomposed 
into the longitudinal and transverse parts, x^il^) = 
X L (q,u)q^q„/q 2 + X T (<7, uj) (S^ - q^/q 2 ) , in which 
uia 2 (iv) belongs to the transverse part. (In the case of 
the rotation, the effect of the wall of the bucket propa- 
gates along the radial direction, which is perpendicular 
to that of particle motion.) Hence, we rewrite eq.(7) as 

pda 1 ( U ) = --J o du u2 _ ua . (8) 

In the normal-fluid phase, because of x T (Qt uj ) = 
X L (q,oo) for small q, x T (q,cj) in eq.(8) is replaced by 
the longitudinal part x L (Qt uj )- Such a modified formula 
with the aid of p = x L (g,0) defines the normal fluid 
conductivity a\ n {uj) (see Appendix. B). 

As T — > T\, owing to Bose statistics, the energy of 
the low- lying transverse excitation increases (see §2.3). 
It destroys the balance between the transverse and the 
longitudinal excitations, hence x L (<Z: UJ )^X T (q, w) 7^ at 
lo — > 0. Instead of a\ n {uj) defined by x L (g,o;), one must 
go back to the original definition of eq.(8), and rewrite 
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Fig. 3. p s {T)/p obtained by eq.(12) using v{T) in Fig.l. 
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ai(u>') = ai„(w') 
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Y\m q ^ Q [x L {q,uj) - x T (q,u)} 



.(9) 



The integral over to in Eq.(9) leads to a sharp peak at 
lo' = 19 ^ with the aid of Hilbert transformation as 



dto 



(10) 



Accordingly, a sharp peak lim 9 _,.o [x L (<2S 0) — 
X T (q, 0)](5(cj) appears in <ti(uj) of eq.(9) as required in 
eq.(3), where lim^o [x L (<7> 0) — X T (<7)0)] corresponds to 
the mesoscopic superfluid density p s (T). Using eq.(10) 
in eq.(9), and using the result in a — uj 1 J^ f a(u)cLj, 
we take out the sharp peak from this integral, and 
obtain the kinematic viscosity v{T) = l/(4er) as 21 ) 



v{T) 



2ir p s {T) 



(11) 



d 2 ojf p ™ 
where v n (— l/[4cr„]) is v of the classical liquid. 



2.2 Comparison with the experiment 

Equation (11) is a general formula depending on no 
specific model of liquids, and we can compare it to v (T) 
shown in Fig.l. 10 ) When we view the fall of v(T) in Fig.l 
as a manifestation of a mesoscopic condensate, we can 
estimate its p s /p using eq.(ll) as follows. We assume v n 
in eq.(ll) to be v at T „=3.7K in Fig. 1, and derive ujj in 
eq.(ll) from the time dependence of the level of a liquid 
h(t) = h(0)exp(—at) as follows. The a of liquid helium 3 
was measured at 1.105 K as 5 x 10~ 4 /s in Fig. 2 of ref.9. 
At 3.7 K, it approximately changes to 3 x 10~ 4 /s with 
a change of p/rj in a oc pgd /(rjL). For liquid heliums 
3 and 4 near 3.7 K, the factor p/n is nearly the same. 
Hence, we estimate that liquid helium 4 at 3.7 K has 
a = 3 x 10~ 4 /s for a typical size of the capillary, and 
LJf = \^3a = 5 x 10~ 4 rad/s (a small value). With the 
above ujf and a typical capillary radius d = 10~ 2 cto, we 




Fig. 4. Deformation of a small rectangular clement in the flow. 
The shear stress F xy is expressed as r]dv x /dy. 



estimate p s (T)/p using v(T) of Fig.l as 



2t: Ps (T) 



= v(T)- 1 - v{T on y 



(12) 



d 2 ujf p 

and obtain Fig. 3. Just above T\, p s (T)/p reaches 1.6 x 
10" 5 . 

Compared with the nondissipative process such as the 
rotation of liquid, the intrinsically dissipative process 
such as the capillary flow depends on some experimental 
parameters such as w/ and d as in eq.(ll). In the rotating 
bucket experiment, the change in the moment of inertia 
I z simply obeys I Z (T) = if [I - p s {T)/p]. In the rota- 
tion experiment by Hess and Fairbank, 22 ) the moment 
of inertia I z above T\ is slightly smaller than the nor- 
mal phase value If. Using these currently available data, 
p,(T x + 0MK)/p = 8 x 10- 5 , and p,(T x + 0.28K)/p = 
3 x 10~ 5 in ref.14. Despite the subtleness inherent in the 
dissipative process, p s (T)/p in Fig. 3 has the same order 
of magnitude as p s (T)/p independently obtained using 
I Z {T). Just above T\, about 1/10 5 of all helium 4 atoms 
participate in the mesoscopic condensate. They are neg- 
ligible in thermodynamic quantities, but manifest them- 
selves in the mechanical response such as shear viscosity. 
In eq.(ll), the effect of p s (T) is amplified by the small 

u f- 

2. 3 Physical picture in coordinate space 

(1) Let us discuss the physical mechanism behind 
eq.(ll). Figure 4 shows a small rectangular element in 
the flow, which initially deforms macroscopically, then 
relaxes to a part of the stationary flow with a relaxation 
time tm- For 77 of the classical liquid, Maxwell obtained a 
simple formula n — Gtm (the Maxwell relation) in anal- 
ogy with solids, where G is the modulus of rigidity. 23 ) 
The macroscopic relaxation in Fig. 4 is a result of the ac- 
cumulation of microscopic relaxations by the excited par- 
ticles. Quantum mechanics states that, in the decay from 
an excited state with an energy E to a ground state with 
Eq, the higher excitation energy E causes a shorter mi- 
croscopic relaxation time r as h/r ~ \E — E \. In eq.(8), 
the left-hand side includes the tm hi o~\ = p/(AGtm), 
whereas the right-hand side includes the many-body ex- 
citation spectrum in x T (<7,w). In this sense, eq.(8) is the 
many-body version of h/r ~ \E — Eq\. (In other words, 
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cq.(8) is another expression of the Maxwell relation.) In 
liquids, slightly different microscopic local structures are 
irregularly arranged , the energies E's of which differ only 
slightly from each other. In the relaxation of one arrange- 
ment to others, the smaller energy difference \E — Eo\ 
leads to a longer r in K/t ~ \E — Eq\. These long r's 
result in a long tm of the macroscopic deformation, thus 
leading to a large shear viscosity r\ — Gtm of the classical 
liquid. 

In liquid helium 4 near 7\, no structural transition is 
observed in coordinate space. Hence, G must be a con- 
stant at the first approximation, and therefore the fall of 
the shear viscosity r\ — Gtm is attributed to the decrease 
in tm- At the microscopic level, this comes from the de- 
crease in r's. In view of K/t ~ \E — Eo\, the decrease 
in r suggests the increase in E. Near T\, it is natural 
to attribute it to the effect of Bose statistics. The rela- 
tionship between the excitation energy and Bose statis- 
tics dates back to Feynman's argument on the scarcity 
of low-energy excitation in liquid helium 4, 24 ' in which 
he explained how Bose statistics affects the many-body 
wave function in configuration space. We will apply his 
explanation to the analysis of shear viscosity. 

(2) Consider the motion of particles in Fig. 4. For ex- 
ample, two particles I and 2, each of which simultane- 
ously starts at (x, y) and (x,y + Ay), move along the 
x-direction to positions V and 2'. The long thin arrows 
represent the displacement of white circles to black ones. 
In the BEC phase, the many-body wave function has 
permutation symmetry everywhere, and all white circles 
are therefore permutable in Fig. 4. At first sight, these dis- 
placements by long arrows seem to be a large change, but 
they are reproduced by a set of small displacements of 
neighboring white circles by short thick arrows as shown 
in Fig. 4. In Bose statistics, owing to permutation sym- 
metry, one cannot distinguish between two types of par- 
ticles after displacement, one moved from near positions 
by a short arrow, and the other moved from distant ini- 
tial positions by a long arrow. Even if the displacement 
made by the long arrows is a large displacement in classi- 
cal statistics, it is only a small displacement by the short 
arrows in Bose statistics . The displacement related to 
shear viscosity is a transverse one. In general, the trans- 
verse displacement does not change particle density on 
a large scale, and therefore, for any given particle after 
displacement, a nearby particle always exists in the ini- 
tial distribution. (On the other hand, the longitudinal 
displacement largely changes the particle density, and 
therefore nearby particles do not always exist, which im- 
plies x L (v,u) - x T (q,u) ^ in eq.(9).) 

Let us consider this situation in 3N-dimensional con- 
figuration space following Feynman. The excited state 
related to shear viscosity, in which particles make small 
displacements, has a wave function that is not far apart 
from the ground-state wave function in configuration 
space. The excited-state wave function must be orthog- 
onal to the ground-state wave function. Since the latter 
has an uniform amplitude, the former must spatially os- 
cillate between the plus and minus values. This means 
that the wave function of the excited state due to slight 
displacements must oscillate within a small distance in 



configuration space. The kinetic energy of the system is 
determined by the 3N-dimensional gradient of the many- 
body wave function, and therefore this steep rise and fall 
of amplitude raises excitation energy. The relaxation of 
such a microscopic state is a rapid process. When these 
processes occur simultaneously, it leads to a rapid relax- 
ation of the macroscopic deformation. This mechanism 
intuitively explains why Bose statistics leads to the de- 
crease in shear viscosity i] — Gtm- 

(3) When the system is at high temperatures, the co- 
herent wave function has a microscopic size. If the long 
arrow in Fig. 4 takes a particle out of such a wave func- 
tion, the particle after displacement cannot be repro- 
duced by the slight displacement. The mechanism be- 
low T\ does not work for such a displacement, and tm 
changes to the ordinary long tm in the classical liquid. 
When the system is in T\ < T < 3.7 K, the size of the co- 
herent wave function is mesoscopic. In a repulsive system 
with high density such as liquid, the large-distance mo- 
tion takes much energy, and therefore in the low-energy 
excitation, particles are likely to stay within the same 
wave function. The energy of such an excitation is low 
compared with that of the excitation due to large dis- 
placements, but owing to Bose statistics, it is not so low 
as that in the classical liquid, and therefore its relax- 
ation to the ground state is relatively a fast process. 12 ) As 
T — >• T\, the number of such fast relaxations gradually 
increases, and the superflow appears within the meso- 
scopic distance, which is the reason for the substantial 
fall of r\ at T A < T < 3.7K. 

3. Microscopic Model 

3. 1 Onset of superflow 

To formulate the mechanism of the fall of v above T\, 
we consider the repulsive Bose system having the follow- 
ing hamiltonian 

H = J2 e( P )^ p %+Uj2J2 K-A'+^p'^ ( u > 
p p,p' q 

(13) 

where $ p is the annihilation operator of a spinless bo- 
son, and begin with the state without the macroscopic 
condensate. The onset of superflow above T\ has a sim- 
ilar mechanism to the onset of the nonclassical moment 
of inertia above T\. In this subsection, we recapitulate 
ref.14 on this point with some modifications. 

Let us discuss the role of particle interaction. The 
susceptibility x^v 

(q, lj) is a correlation function of cur- 
rents j^q,T) = E P ,n(P + ?/2))1>JV 9 e ^" T - From 
Xiiv = X L (l^u/q 2 + x T {^ v - q^q v /q 2 ), we will extract 
the term proportional to q^q v , and define x^^u;) = 
(x L — X T )QhQv/q 2 f° r Ps(T). In the ideal Bose system, 
we obtain 



1 M e (p)) - fs{e{p + q)) 



4 V 



E 



u + e(p) - e(p + q) 



(14) 

where /s(e(p)) is the Bose distribution. If bosons would 
form the macroscopic condensate, /s(e(p)) in eq.(14) is a 
macroscopic number for p = and nearly zero for p ^ 0. 
Thus, in the sum over p on the right-hand side of eq.(14), 
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Fig. 5. When p = p' and p + q = p' + q' in (a), an exchange of 
particles between large and small bubbles yield (b). 



T\ 3 To 




Fig. 6. The temperature dependence of fi(T) — £(f/) depicted by 
a solid curve, and the condition of eq.(20) by a one-point-dottcd 
one. T on corresponds to 3,7K in Fig.l. 



only two terms corresponding to p = and p = —q 
remain, with a result of x^g, 0) = p s (T)q fl q„/q 2 . When 
bosons form no condensate, however, the sum over p in 
eq.(14) is carried out by replacing it with an integral, 
and one notices that q~ 2 dependence disappears, hence 
p s (T) = 0. This means that when examining the system 
near T\, under the particle interaction Hj is 

needed. We obtain a perturbation expansion of X^v{q, w ) 
with respect to U as 



(G\T T jn(x,T)j„(0, 0)\G) 



(15) 



(0\T T j^(x,T)j v (0,0)exp 


- / M ; (r) 
Jo 


|0> 


(0 \exp 


- / drH^r) 
Jo 


10) 



In eq.(8) and its modified form of <7i„(w) using 
X L {q,u), the characteristic excitations of a liquid is in- 
cluded in x L {q,u). X^(g,w) = (x L - X T )%1v/q 2 repre- 
sents the change of x L induced by Bose statistics. When 
Xnu{q,u) first obtains a nonzero value at T oni one can 
represent it by the simplest process in eq.(15). Figure. 
5(a) shows j^j v as a large bubble, and the effect of Hi 
as a small bubble with a dotted line U. As T — > T\ in 
the normal phase, Bose statistics forces particles in the 
large and small bubbles to strictly obey the permuta- 
tion symmetry. In Fig. 5(a), when one of two momenta 
are equal (p = p') for two bubbles, and when the other 
two momenta are also equal (p + q = p' + q'), a graph 
made by exchanging these particles must be included in 
the perturbation expansion. Cutting the line at the point 
denoted by opposite arrows in Fig. 5(a), and reconnect- 
ing it to the line with the same momentum of the other 
bubble yields two bubbles in Fig. 5(b) with the common 
momenta p and p + q connected by the repulsive interac- 
tion U. As the mesoscopic condensate grows, such an ex- 
change occurs many times, thus leading to a large chain 
of bubbles. Furthermore, among various momenta, the 
process including particles with p = 0oip + q = grows 
to play a dominant role. As a result, we obtain 



(16) 



l 



4 V 



Fp(q,iuj) 



Fp(q,-iu}) 



l-UF^{q,iuj) l-UF fj (q,-iuj) 



where 

F p (q,iu)= (17) 
(expQ3[g - p}) - I)' 1 - (cxp(/3[e(q) + E- /i]) - l)' 1 
-ico + e(q) 

(E is the self energy and [i is the chemical potential.) 
Since our interest is x L ~ X T near w = 0, we expand 
eq.(16) with respect to uj 



x^(q,u) 



q^q u 1 2Fp(q,0) 
4 Vl-UFp(q,0) 



e(q)(l-UF (q,O)) 



+ 



(18) 



Fp(q, 0) in eq.(17) is a positive monotonically decreas- 
ing function of q 2 . Since our interest is the macroscopic 
response of the system, let us expand it with respect to 



q 2 as 







4 S inh^(«^ 



1 - 



= a — bq 



1 

^anh(M^h2m 



(19) 



The denominator on the right-hand side of eq.(18) has 
the form 1 - UF p (q, 0) = (1 - Ua) + Ubq 2 H . As T ap- 
proaches T\ (fj,(T)—T, — > 0), the first term a in the expan- 
sion of Fp(q, 0) gradually increases. Hence, 1 — Ua grad- 
ually decreases, and finally satisfies 1 — Ua = 0. At this 
temperature, eq.(18) becomes proportional to q^q^/q 2 as 
V r_1 [a/ (2Ub)]q fl q l/ /q 2 , hence giving its coefficient x L ^X T 
a nonzero value. This is a sign of the onset of superflow 
by the mesoscopic condensate . 14 ) 
The condition 1 — Ua = yields 



Uf3 = 4sinh" 



(20) 



Figure. 6 schematically shows the condition of eq.(20) 
by a one-point-dotted curve, and /x(T) — S of the system 
by a thick curve. (The simplest approximation for /u(T) 
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of liquid helium 4 is to use /u(T) of the ideal Bose gas 
with T c replaced with T\, 



H(T) - E([7) 



.93/2(1) 
20? 



rp \ 3/2 



(21) 

which approximation dates back to London's paper.) 
Since /i(T) — E reaches zero at the finite temperature 
T\, it approaches zero more rapidly than the one-point- 
dottcd curve by eq.(20). Hence, above T\, there is al- 
ways a temperature in which eq.(20) is satisfied. This is 
the reason why v(T) in eq.(ll) begins to fall above T\ 
as in Fig.l. We call this temperature the onset tempera- 
ture T on . Using the definition of a and b in eq.(19) and 
Ua = 1 in V~ 1 [a/(2Ub)]q a q v /q 2 , the mesoscopic super- 
fluid density p s (T) = lim^o [x L ~ X T ] at T on has the 
form 

I m 

P ' s (T) = V sinh|/3[//(T)-E]| ' (22) 

Since the sizes of the condensate and capillary are not 
so different, we do not take the V — > 00 limit in eq.(22). 
Hence even when /x— E 7^ 0, p s (T) is finite. Using eq.(22) 
in eq.(ll), we can describe the fall of v(T) in Fig.l. Equa- 
tion (22) serves as an interpolation formula of p s (T) in 
T\ < T < 3.7 K, and jumps to a macroscopic value at 
T x . 

3.2 Conductivity spectrum 

Let us obtain the expression of cr 1 (w) in eq.(9). In 
Fig. 2(b), the emergence of the peak at u> = is accom- 
panied by a gradual change of the continuum part. At 
T = T on (Ua = 1), we rewrite eq.(18) as 



where 



q 2 2Vb 



1-^ 

w5 



4 ^ T )\anh^^ 



(23) 



(24) 



(q e is the wave number representing the size of the 
capillary.) We want a simple formula of Xjw(9> w ) that 
agrees with eq.(23) at u — > and shows a reason- 
able behavior at u — > 00. At T „, the simplest form of 
= (x L — X T )q^qv/ q 2 satisfying these conditions is 
given by 

.2 1 

(25) 



~ / \ q»qv a 

Xa"A<2S w ) - — 



q 2 2Vb 



1 " 



Equation. (25) represents the dynamical balance between 
the longitudinal and transverse responses. Using the def- 
initions of a and b in eq.(19) and Ua — 1, we define 
the generalized form of the mesoscopic superfluid den- 
sity Ps{T,u>) = \im q ^ [x L (q,uj) - x T (<7,w)] as 

MT, U ) = ^ —^= =, , * ■ (26) 



Usinh|/3[/i(T)-E]| 



1 + 



Using eq.(26) as \ L ~X T in °q-(9), taking the finite part 
of the integral as 



duj- 



1 



1 



— 7T 



(27) 



2w d 



and eq.(24) for u;<j, we obtain the real part of conduc- 
tivity as 



ci (w) = cti„(w) + 



2U- 



pird 2 sinh|/3[/x(T) - E]| 2 



■(5(cj) (28) 



2U" 1 / k R T 



pud 2 \e(q e ) 2 



cosh' 



^[Mr)-E]| 



1 + 



7? 



<j\ (w) must satisfy the sum rule eq. (5) whether it is in the 
normal or superfluid phases. (In ai(uj)duj of the sum 
rule, the second term in the bracket of the right-hand side 
of eq.(28) yields a term proportional to 1/ sinh \f3[p(T) — 
E]| with the aid of J °° duj /(l + uj 2 /uj^) = irujd/2 and 
eq.(24). The first and second terms in the bracket of 
eq.(28) cancel out each other in J °° ai(uj)dio. 21 ^) 

In eq.(28), <7i„(w) is given by the real part of eq.(B5) 
in Appendix. B. Hence, the conductivity of the normal- 
fluid part a n (uj) is given by 

/ \ 

1 



to- 

2 
nd 2 



V 



k B T 

<q e ? 



./ (1 (id[l+, lv f_ 

v- 1 



(29) 



cosh' 



)(- 








<4) 



(Jo is the zeroth order Bessel function, and n = p/m). 
The temperature dependence of cr n (uj) comes from the 
second term on the right-hand side of eq.(29). As fi(T) — 
E -5- 0, 1/cosh 2 \(3[p(T) - E]/2| and increase, hence, 
Jo° o- n (uj)duj decreases. On the other hand, the conduc- 
tivity of the superfluid part a s (oj) satisfies 



f 

Jo 



<j s (u))du) 



V~ 



(30) 



2nd 2 sinh|/3[/j(T) - E]|' 

As p{T) — E — > 0, the sharp peak at w = increases 
owing to l/sinh|/3[/x(T) - E]|. Equation (29) and (30) 
are the simplest formulae describing the change in a(co) 
at T x < T < 3.7K in Fig.2. 

3. 3 Stability of superflow 

Equation. (26) implies that, for the objects oscillating 
at the frequency above wa, a considerable part of the su- 
perfluid behaves as a normal flow. Here, we call u>d in 
eq.(24) the damping angular frequency . For the exper- 
iment using the vertically standing capillary in gravity, 
we measure the integrated response of the system from 
u> = to ujf. As long as the properties at approximately 
lo = determine the observed quantities, this feature 
above ujd does not seriously affect the result. On the other 
hand, there is another type of experiment in which the 
superfluid behavior is measured only at certain frequen- 
cies, such as torsional oscillation or ultrasound. In this 
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method, there is a case in which the above feature gives 
rise to a serious problem (see §4). 

The stability of the superflow by the mesoscopic con- 
densate depends on temperature. As the Bose statisti- 
cal coherence grows in eq.(24) (p(T) — E — > 0), the 
superflow becomes robust to the oscillating probe in 
eq.(26) (ujd — > oo). By gradually changing temperature at 
T > T\, we can realize different size distributions of the 
mesoscopic condensate in the capillary flow. The quan- 
tity reflecting the size of the condensate is the chemical 
potential p,(T). Here, to characterize the size distribution 
of the condensate at T, we use the number of particles 
(s) satisfying — E)(s) = —1 (see Appendix. C). For a 
small fi — E, eq.(24) is rewritten as 

-=W<»>- (3I) 

As (s) grows, ujd increases, hence, as the size of the meso- 
scopic condensate grows, the superflow by the mesoscopic 
condensate gradually becomes robust to the oscillating 
probe. 

The stability of superflow is also related to the re- 
pulsive interaction U. In the famous argument by Lan- 
dau on the critical velocity of fully developed BEC, the 
change of one-particle spectrum from p 2 /(2m) to v s p 
induced by the repulsive interaction in the Bogoliubov 
spectrum plays a crucial role. 2 ) Although Landau's argu- 
ment deals with the magnitude of critical velocity, eq.(26) 
deals with the upper limit of frequency at which flow 
velocity changes. When U = in cq.(18), p s (T,uj) is 
written in the form of eq.(26), with sinh |/3[/x(T) — E]| 
replaced by sinh 2 (|/3[/z(T) - E]|/2). In this case, lo a 
in cq.(24) is simply e(q e ) 2 . Hence, with increasing lj, 
Ps(T,lo) vanishes far more rapidly than that in the case 
of U 7^ 0. Physically, the repulsive interaction U prevents 
the drop of particles from the superflow, thus stabilizing 
it. In eq.(24), the repulsive interaction gives the factor 
I/tanh 2 (| / 3[ / u(T)-E]|/2) to uj 2 . Hence, p s {T,u) does not 
easily vanish for a large w in eq.(26). The repulsive inter- 
action plays the significant role not only in the emergence 
of superflow by mesoscopic condensate, but also in its sta- 
bilization. In this sense, eqs.(24) and (26) are extensions 
of Landau's argument to the underdeveloped BEC. 

4. Implication for Superflow in Porous Media 

In the capillary flow above T\, the mesoscopic con- 
densate spontaneously emerges in free space, whereas in 
porous media below T c , it is forced to emerge in a re- 
stricted space. 25 ) Despite this difference, when we focus 
on a single flow through a pore, the result of this paper 
gives us a simple criterion for the stability of superflow in 
porous media. In porous media, the superflow is normally 
measured by the method of torsional oscillation or ultra- 
sound. Just below T c , (s) in eq.(31) and ujd in eq.(26) 
are still small. When the oscillator at u> > u!d induces 
the oscillation of superfluid, [p s (T, 0) — p s (T,u)}/ p s (T, 0) 
of the superfluid behaves as a normal- fluid part. Hence, 
w 2 /(w 2 +W 2 .) of superfluid is locked to the substrate , and 
p s (T) is underestimated. With decreasing temperature, 
this underestimation of p s (T) is improved by an increase 
in u>d- 2 ^ Alternatively, when we change the frequency 



at a given temperature, the oscillation methods using a 
higher frequency will detect a larger discrepancy between 
T c in the mechanical response and T c in the thermody- 
namic one. 

On the other hand, in the oscillation experiments on 
the superfluid, another type of upper limit uj p is known, 
which is set for the normal-fluid part to be locked to the 
substrate. (The viscous penetration depth S = vj (2uj) 
must be greater than the channel size. Since the relevant 
size of the channel is about 10 -5 cra and v ~ I0~ 4 cm 2 /s, 
f p = u> p /2tt is roughly 10 6 Hz.) In normal situations, 
Lo p > ujd is usually satisfied. At uj > uj Pi neither the 
normal-fluid nor superfluid parts is locked to the sub- 
strate, hence, p s (T) is overestimated. At ujd > ui, the 
normal and superfluid parts are correctly discriminated. 
A delicate problem arises at uj p > uj > u>d, especially 
when a frequency comparable to f p = uj p /2tt is used 
in ultrasound. Although ui p does not depend on tem- 
perature, LUd increases with decreasing temperature, and 
there may be a case of ujd > uj p at very low temperatures. 
When we represent ujd by its maximum ujd(T = 0), we 
can classify the situation into two cases: uj p > uJd(T = 0) 
and uJd(T = 0) > cj p , depending on the type of porous 
medium. The interpretation of the experimental data ob- 
tained in such a situation will be complicated, because 
of the existence of two different kinds of upper limit in 
frequency. 

Recently, a comparison between the measurement us- 
ing a torsional oscillator (2140 Hz) and by ultrasound (10 
MHz) was made in liquid helium 4 confined in a porous 
substrate, hectorite, and Gclsil. 27 ) The decoupling of the 
superfluid from the oscillating substrate results in a de- 
crease in the moment of inertia AI Z , and in the increase 
in the ultrasound velocity Av s . Different situations will 
be realized according to different frequencies and sub- 
strates. The comparison of experimental data between 
different frequencies will give us a clue to the dynamics 
of superflow in porous media. 

5. Discussion 

5.1 Comparison between superflow in the dissipative 
process and that in the nondissipative one 
On the onset mechanism of superflow, it makes a dif- 
ference whether superflow appears in the dissipative or 
nondissipative process. Comparing the moment of in- 
ertia I Z {T) = I z (1 — p s (T)/p) for the non-dissipative 
process 14 ' with the kinematical viscosity v{T) of eq.(ll) 
for the dissipative one, we note the following features of 
u(T). 

(1) In I Z {T) = if (I - p s {T)/p), p s (T) appears only in 
the coefficient of the linear term of if. In v(T) of eq.(ll), 
the effect of Bose statistics appears in all coefficients of 
higher-order terms of v n except for the first-order one. 
This feature does not depend on the particular model 
of a liquid, but on the general argument. (On the other 
hand, the microscopic derivation of v n depends on the 
model of a liquid.) 

(2) In I Z (T) = If {I - p s {T)/p), the change in 
X L {q,<jj) — x T (g,w) directly affects I z without being en- 
hanced, and therefore the observed effect of the finite 
Ps(T) above T\ is very small. In v{T) of cq.(II), because 
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of the sharp peak in the dispersion integral in eq.(9), 
the small change in p s {T) is strongly enhanced to the 
observable change in v(T). 

(3) The existence of 1/cP before p s (T)/p in eq.(ll) 
indicates that the narrower capillary shows clearer evi- 
dence of frictionless flow. Similarly, the existence of 1/uJf 
indicates that the choice of experimental procedure, such 
as the method of applying the pressure between two ends 
of the capillary, affects the temperature dependence of 
the v{T) of capillary flow. This means that superflow ap- 
pearing in the dissipative process depends on more vari- 
ables than that in the non- dissipative one, which is in 
accordance with the general feature of nonequilibrium 
phenomena. 

5.2 Comparison with thermal conductivity 

In liquid helium 4 near T\, we know a marked change 
in another type of conductivity, the anomalous thermal 
conductivity. Under a given temperature gradient VT, 
the heat flow Q satisfies Q = — kVT, where n is the 
coefficient of thermal conductivity. In the critical region 
above T\ (\T/T X - 1| < 1CT 3 ), the rapid rise of k is 
observed, and finally at T — T\, n abruptly jumps to 
infinity. For T\ < T < 3.7 K, however, two kinds of con- 
ductivity behave differently. In shear viscosity, the cor- 
responding conductivity a(T) = 1/(4^) shows a gradual 
rise, whereas n{T) shows no such rise. This difference is 
expected for the following reason. The k in thermal con- 
ductivity is expressed in terms of the correlation func- 
tion with a similar structure to the r\ in shear viscosity, 
but k and r\ are qualitatively different. Although shear 
viscosity is associated with the transport of momentum 
(a vector), thermal conductivity is associated with that 
of energy (a scalar). For the vector field (velocity field), 
the direction of vectors has a rich variety in its spatial 
distribution. Among various flow-velocity fields, we can 
regard the Couette flow rotating like a rigid body in a 
bucket as the nondissipative counterpart to the capillary 
flow. On the other hand, for the scalar field (temperature 
field), the variety of possible spatial distributions is far 
limited. The flow of heat energy is always a dissipative 
phenomenon, and there is no nondissipative counterpart. 
Hence, the formulation using the Kramers-Kronig rela- 
tion in §2 cannot be applied to the onset of the anoma- 
lous thermal conductivity, and no mechanism that am- 
plifies the small p s {T) is expected. This formal differ- 
ence between shear viscosity and thermal conductivity 
is consistent with the experimental difference between 
a(T) = l/(4z/) and k(T) at T A < T < 3.7K. 

5.3 Comparison with Fermi liquids 

The fall of the shear viscosity in liquid helium 3 at 
T c is a phenomenon parallel to that in liquid helium 4. 
The formalism in §2 is applicable to liquid helium 3 as 
well. For the behavior above T c , however, there is a strik- 
ing difference between liquid helium 3 and 4. The phe- 
nomenon occurring in fermions in the vicinity of T c is not 
the gradual growth of the coherent wave function, but the 
formation of Cooper pairs by two fermions. (This differ- 
ence evidently appears in the temperature dependence of 
specific heat: The C(T) of liquid helium 3 shows a sharp 



peak at T c without a symptom of its rise above T c .) Once 
Cooper pairs formed, they are composite bosons with a 
high density at low temperatures, and immediately jump 
to the superfluid state. Hence, the shear viscosity of liq- 
uid helium 3 shows an abrupt drop at T c without a grad- 
ual fall above T c , like the dotted line in Fig.l. 

5-4 Future problems 

As T approaches T\, eq.(22) becomes insufficient to 
use as an interpolation formula of p s (T) , because it pre- 
dicts a far smaller p s /p than that in Fig. 3 when eq.(21) is 
used. This implies that a considerable number of helium 4 
atoms with p^O are dragged into the superflow. Owing 
to the repulsive interaction, particles are likely to spread 
uniformly in coordinate space. By this feature, particles 
with p are forced to behave similarly to particles 
with p = 0. // they move differently from the superflow, 
particle density becomes locally high, thus increasing in- 
teraction energy. This is the reason why a considerable 
number of particles with p ^ participate in the super- 
flow even when it is above T\ . This tendency will become 
apparent as T — > T\. Hence, the observed p s (T) may be 
the sum of all p s (T) 's over different momenta. To explain 
Fig. 3 quantitatively, a more realistic model of liquids is 
needed. 

Among all nonequilibrium phenomena in liquids, the 
capillary flow of liquid helium 4 takes a unique position. 
When increasing temperature from below to above T\, 
one can see a primitive form of the shear viscosity of liq- 
uids in the gradual rise in v(T), which arises from the 
superfluid in such a way that it does not occur in other 
systems. The emergence of slightly different local struc- 
tures leads to a long relaxation time tm, and hence to a 
finite shear viscosity of liquids. Its study has a possibility 
of shedding new light on the liquid theory. 

In this paper, we used the Poiseuille solution of clas- 
sical fluid dynamics, and examined the quantum effect 
appearing in the coefficient of shear viscosity. As an al- 
ternative approach, we know the quantum hydrodynam- 
ics by Landau 2 ). 28 ' If we consider the shear viscosity at 
T\ < T < 3.7K using this approach, it will correspond 
to the semiclassical solution in quantum hydrodynamics. 
This will give us interesting problems from practical and 
formal viewpoints. 

Appendix A: h(t) — h(0) exp(— at) 

To measure r\, let us set the capillary to stand vertically 
in gravity g, the upper end of which is connected to a 
reservoir, while the lower end is open to the helium bath 
(Fig.l in ref. 8). The mass of liquid passing through the 
capillary per unit time is Q — 2irp J Q rv x (r)dr, which is 
given by eq.(l) as 



h(t) is the difference between the level inside the reser- 
voir and the level in the bath. In the reservoir with a 
radius R, a mass of liquid pnR 2 h(t) flows out through 
the capillary at a rate of Q. In the vertically standing 
capillary, AP(t) between the above two levels is pgh(t) 
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in eq.(Al). h(t) decreases according to 



, dh 



dt 



8?7 L 



(A-2) 



thereby leading to h(t) = h(0) cxp(— at) with a = 
pgdy(8 V LR 2 ). 

Appendix B: Conductivity Spectrum cr(u)) 

The Stokes equation under the oscillating pressure gra- 
dient AP cxp(iujt) I L is written in the cylindrical polar 
coordinate as 

dv ( d 

dt = 

Velocity has the form 



d \ APexp(iwi) 
dr 2 rdr J pL 



(B-l) 



v(r,t) 



AP exp(iut) 



iujpL 

under the boundary condition v(d,t) 
Av(r,t) satisfies 



+ Av(r,t), 



(B-2) 



dAv(r,t) 
dt 



d 



d 



dr 2 rdr 



0. In eq.(Bl), 
Av(r,t), (B-3) 



which solution is written in terms of the Bessel function 
Ja(iXr) with A = (I + i)\Ju/ (2v). Hence, 



v(r,t) = 



APexp(iwi) 
iujpL 



1 - 



Jo (iXd) 



(B-4) 



Using v(r = 0,t) of eq.(B4), the conductivity spectrum 
a(u) satisfying pv(0,t) = cr(uj)d 2 AP exp(iwi) / L (eq.(2)) 
is given by 

/ \ 



1 



icud 2 



1 



V 



Jo 



id[i+ t y— 



(B-5) 



The real part of eq.(B5) gives a curve of cr(u>) in Fig. 2(a). 
(Re <t(0) in eq.(B5) agrees with p/{Ar)).) It also deter- 
mines the conserved quantity ir^ 1 / <r(uj)duj = f(d) cx 
d~ 2 in eq.(5). 

Appendix C: Size of the Condensate 

The grand partition function of the ideal Bose sys- 
tem Z (/j) = Y\(l — cxp[/3(e p — p)])~ x was rewritten by 
Feynman and Matsubara in terms of the size distribution 
of the coherent wave function. 29 ) For the repulsive Bose 
system with p,(T) — E, it is given by 



a rough estimate of the typical volume size of the con- 
densate. 
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